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Abstract 



We study the sample-size dependence of the ground-state energy in a 
one-dimensional localization problem, based on a supersymmetric quantum 
mechanical Hamiltonian with random Gaussian potential. We determine, in 
the form of bounds, the precise form of this dependence and show that the 
disorder-average ground-state energy decreases with an increase of the size 
R of the sample as a stretched-exponential function, exp(— R z ), where the 
characteristic exponent z depends merely on the nature of correlations in the 
random potential. In the particular case where the potential is distributed as 
a Gaussian white noise we prove that z = 1/3. We also predict the value of z 
in the general case of Gaussian random potentials with correlations. 
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I. INTRODUCTION 



The one-dimensional Schrodinger Hamiltonian 

H = - ' + ^(,) + ^ (1) 
dx z dx 

arises in such diverse areas of quantum mechanics as studies of solitons in conjugated poly- 
mers (polyacetylene) fl]||, in which 4>(x) describes the dimerization pattern of the carbon- 
hydrogen sequence, or of electrons dynamics in two-dimensional systems subjected to random 
magnetic fields ||. It also represents a celebrated toy model of supersymmetric quantum 
mechanics introduced by Witten [|J. 

In case of non-random potentials <f>{x) the properties of the Hamiltonian in Eq.(l) have 
been thoroughly studied within the last two decades and several interesting results have been 
obtained. In particular, exact solutions of the Schrodinger equation have been derived for a 
class of the so-called shape invariant potentials || . Besides, such form of the Hamiltonian 
inspired a new method of semi-classical quantization [§,0. 

More recent studies of H have focused on situations in which <p(x) is random |3|,|8HTll. 
Here Eq.(l) describes a localization problem in which, remarkably, the average density of 
states, the localization length and the Lyapunov exponent can be computed exactly |^-|10|. 
In fact, this problem represents an additional to the two known by now |12|-|l4| examples of 
localization problems which can be solved exactly in continuum. 

In the present paper we study a new aspect of the one-dimensional localization problem 
associated with the Hamiltonian in Eq.(l) in which the potential <p(x) is a random function 
of the space variable x. We focus on the physically interesting situation where H is defined 
on a finite interval [-R, R] and analyse the sample-size dependence of the ground-state 
energy E (R, {(f)}), which may be expressed in terms of non-local exponential functionals of 
random potential <p(x). We consider mostly throughout the paper the case where <j>(x) is 
Gaussian, 5-correlated white noise with the moments 

< > = (2) 
1 



and 



< (f)(x)(j)(x') > = cr6(x - x'), (3) 

where the brackets denote the average with respect to realizations of 4>(x). The relevance 
of short-range correlations in the distribution of 4>(x) and their influence on the sample-size 
dependence are also succinctly addressed and several predictions are made. We find that for 
such potentials the disorder-averaged ground state (DAGS) energy 

E {R) =<E (R,{(P})> 

decreases as the size R grows as a stretched-exponential function, E (R) oc exp(— R z ), 
where the characteristic exponent z depends only on the correlation properties of the random 
potential. In the particular case where fluctuations in 0(x) are delta-correlated, as in Eq.(3), 
we prove that z = 1/3. The value of the exponent z is also predicted, on heuristic grounds, 
for the general case of random Gaussian potentials with correlations. We show that such a 
stretched-exponential dependence on R stems from atypical realizations of <p(x), which are 
reminiscent of the representative trajectories supporting long-time anomalous decay laws of 
the survival probability for diffusion in the presence of randomly placed traps |T5| , |TE| or of 
the Lifschitz singularities [jnj] in the low-energy spectrum of an electron in the prese nee of 



randomly dispersed scatterers. Our results are presented in form of lower and upper bounds 
on E (R), which show the same dependence on R and thus define this dependence exactly. 
The method of the derivation of bounds, which we invoke here, has been previously discussed 



in |T8| , |T9l and is based on the statistics of extremes of random potential (j)(x). 

The paper is structured as follows: In Section 2 we present a simple derivation of the 
ground-state energy E (R,{(p}) in case of deterministic potentials 4>(x). We recover in a 
straightforward way an explicit formula for E (R, {(f)}), which coincides when the appropriate 
notations are introduced with the one found in ||20|| , where the influence of the finite-size 
effects on the ground-state energy in conventional quantum mechanics has been examined. 
We then present some arguments showing that this formula is still valid in the case of a 



random potential <f>(x). In Section 3 we estimate the sample-size dependence of the ground- 
state energy considering only typical realizations of random potential (f>(x). Employing the 
standard Jensen inequality we show then that such an estimate constitutes a lower bound 
on the average ground-state energy. Additionally, we illustrate that such an estimate allows 
to recover the correct low-energy behavior of the integrated density of states of Hamiltonian 
(1). Further on, in Section 4, we devise a more accurate approach and derive a lower bound 
(subsection A) and an upper bound (subsection B) on E (R), which determine exactly its 
behavior in the large- R limit. In subsection C we address the question of the DAGS behavior 
in situations, in which the fluctuations of the random potential (f>(x) are correlated, and 
also discuss some similar features between the realizations of disorder supporting large-i? 
behavior of the DAGS and the realizations of random walks which support anomalous long- 
time tails of the survival probability for the diffusion in the presence of traps. Finally, in 
Section 5 we conclude with a summary of our results. 

II. CALCULATION OF THE ENERGY SHIFT IN A FINITE SAMPLE. 

The structure of the Hamiltonian in Eq.(l) is such that, for an arbitrary function <j)(x), 
two independent solutions of the Schrodinger equation 



H^ 2 \x) 



0, (4) 



may be explicitly expressed as functionals of (j)(x) : 




and 
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A. Deterministic potentials 4>{x). 



In this subsection we consider the case of deterministic potentials (f>(x), which decay fast 
enough when x — > ±00 to make ip^ normalizable on the whole axis, i.e. 

/ dx[(po~\x)] 2 < 00 
J — 00 

Therefore, (p^ is a zero-energy wave-function of the Hamiltonian H defined on the whole 
axis. 

Consider now how the situation will be changed if one assumes that the Schrodinger 
equation with the Hamiltonian in Eq.(l) is defined not on the whole X-axis, but only on 
a finite interval [—R,R]. The new ground-state wave-function ^o(x) on this finite interval 
satisfies the Schrodinger equation 

h ^ + <f> 2 (x) + ^] *«,(*) = MR, W) *o(*), (7) 

with an a priori unknown energy E (R, {0}) that will depend on the explicit form of the 
potential (f>(x). Eq.(7) has to be supplemented by the following Dirichlet boundary conditions 
at the ends of the interval, i.e. at points x = —R and x = R, 

^ (x = -R) = ^ (x = R) = (8) 

Our goal will be to estimate the energy shift E (R, {0}) of the ground state caused by the 
introduction of boundary conditions imposed on a finite interval. 

Multiplying both sides of Eq.(7) by (p^\x) and integrating from — R to R, we get the 
following identity 

£ o (iU0}) [ R dx^\x)Mx) = 

J-R 

= J R R dx^\x) h ^ + 2 (*) + ^] (9) 
Integrating by parts the kinetic term on the right-hand-side of Eq.(9) yields 
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E (R, {((>}) [ R dx<pP(x) * (x) = 

J-R 

= &{-K> ^ - rf'(fl) ^ (10) 

ax |x=-_r ax \ X =R 

To proceed further on, we assume that R is sufficiently large and estimate the behavior 
of the derivative of &o(x) i n this limit. Since E (R, {0}) vanishes when R — > oo, one expects 
that the ground-state wave-function on the finite (large) interval may be well approximated 
by a suitable linear combination of the two independent zero-energy solutions, ip^ (x) and 

(2) 

ip (x). In the vicinity of x = R one has 

¥„(*) « rf^) [1 + a [-1^—] (11) 

•/o [<^ ; (x)] 2 

where a is determined through the boundary condition at x = R, which gives 

7o [rfV)] 2 

and consequently, we find that the derivative of &q(x) obeys 

d^o(x) a m f R dx 



Similarly, for the derivative of the ground-state wave-function in the vicinity of x = —R 
one finds 

dx „.-r l ™ Jo [v<o\x)Y 

Now, combining Eqs.(10) to (14) we have 

and, finally, replacing in the integral on the left-hand-side of Eq.(15) the function ^o(^) by 
ip^\x) (which results in exponentially small with R errors pPf ) we find the following explicit 
expression for the ground-state energy on a finite interval 



+ tf[<P { o\x)]- 2 dx J R M\x'W dx> (16) 
Equation (16) reproduces the result of [pOR , derived in terms of a different approach, 
when the notation V(x) = <fi 2 (x) + dcj){x)/dx is introduced. We note also parenthetically 
that expressions of quite identical structure were recently obtained for the diffusion constant 
for random motion in an external periodic potential pT|- f24[| . Therefore, the results which 
will be obtained in the following also apply to this problem, provided that the potential is 
defined as in Eqs.(2) and (3). 



B. Random potentials 4>(x). 

We now have to explain why Eq.(16) is still valid in the case of random potentials <fi{x). 
We first mention that the analysis of the previous subsection is based on the assumption 
that the wave-functions are normalizable on the whole line. In case of random potentials 
4>(x), as defined in Eqs.(2) and (3), the integral J R dx<j)(x) may show an unbounded growth 
when R — > oo and thus the wave-functions are not normalizable. Therefore, it is not possible 
to use directly the results of the previous subsection. One may, however, notice that due 
to the presence of disorder the zero-energy solutions are localized. This observation allows 
us to get rid off the "bad" behavior at infinity. For a given realization of random potential 
4>(x) one chooses R in such a way that 

(j){-R) > and <p(R) < (17) 

and then define an auxiliary configuration 4>(x) (see Figs. la and lb) such that 

<p(x) = (f>(x) for - R < x < R, (18) 

0( x ) = (f>(-R) for x < - R, (19) 

and 

4>(x) = <j>{R) for x > R, (20) 
6 



The corresponding wave-function ^o(%) = exp(/o 4>(x')dx') is therefore an exact zero mode 
of the Hamiltonian 

~ d 2 ~ 9 . . dd)(x) , , 

on the whole line. When i? is sufficiently large, it follows from our previous discussion that 
there exists a ground-state wave-function of the Hamiltonian H on the interval [—R,R], 
whose ground-state energy is given by Eq.(16). Since this state has its support on [—R,R] 
the functions ^q{x) and (po(x) coincide, and thus it is also a quasi zero mode of H with the 
ground-state energy defined by Eq.(16). 

III. TYPICAL REALIZATIONS OF DISORDER AND THE JENSEN 

INEQUALITY. 

We start our analysis of the behavior of the disorder- average ground-state energy in 
Eq.(16) considering first the typical realizations of the random potential <j>(x). To do this, 
we rewrite the ground-state energy for a given sample with a particular realization of the 
random potential <f>(x), Eq.(16), in terms of the following exponential functionals of the 
potential <f>(x): 

T±(z,y) = £ dx [ip { o ] (x)] ±2 = £ dx ex.p(±2 £ <f)(x')dx'), (22) 

which gives 

= 7^^(74^) + 7jhv) (23) 

We notice now that the function W(x) = f$ dx<j)(x), which appears in the definition 
of r_(0, R), Eq.(22), for the potentials as in Eqs.(2) and (3) is simply a trajectory of a 
symmetric random walk. Consequently, for typical realizations of the random potential 
4>(x) the value of W(R) should be of order (aR) 1 ^ 2 ; hence r_(0, R) should grow typically 
as exp(2(cr J R) 1 / 2 ) and then Eq.(23) entails E (R) oc exp(-4( ( x J R) 1 / 2 ). Therefore, for most 
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realizations of <f)(x) one may expect that Eq(R, {<p}) vanishes with an increase of R as 

E (R, {<P})t yP ~ exp(-a 1 / 2 J R 1 / 2 ) (24) 

This typical behavior may be used to evaluate correctly the low-energy asymptotic behav- 
ior of the density of states of the Hamiltonian in Eq.(l). Eq.(24) means that a wave-function 
of low energy E has typically a spatial extension 2R such that 

\n 2 (E) , . 
R oc — — (25) 
a 

Therefore, the number of such states per unit length behaves typically as 



1 a 
2R K W\E)' 



N(E) oc — oc i 2 — , (26) 



It is now interesting to compare Eq.(26) with the exact result §|: 

N ^ = % nwhwr (27) 

where z = \J~Ejo and Jo, Nq are Bessel functions. In the limit E — > one has from Eq.(27) 
that N(E) ~ 2aln~ 2 (E/4:<j 2 ), i.e. the behavior which is quite consistent with our estimate 
in Eq.(26). Therefore, Eq.(26) shows that anomalous singular behavior of the integrated 
density of states is supported by typical realizations of disorder and thus is quite distinct 
from the Lifschitz singularity flI7| , which is most often encountered in disordered quantum 
mechanical systems. 

We now show that the estimate based on typical realizations of the disorder represents 
actually a lower bound on the DAGS energy. To show it explicitly we invoke the stan- 
dard Jensen inequality between the average of the exponential of some function F and the 
exponent of the averaged value of F, 

< exp[-F] > > exp[— < F >] (28) 

Now, choosing 

F = -ln£ o (#,{0}) (29) 



we find, taking into account Eq.(23), that Eq(R) is bounded from below by 
Eq(R) > exp{<ln(r„(-R,R)) > - < ln{r+{-R, R)) > 

- < ln(T-(-R,Q)) > - <ln(r-(0,R)) >] (30) 

One may readily notice that for any random function 4>{x) of zero mean, not all terms in 
the exponent on the right-hand-side of Eq.(30) are to be calculated independently; obviously, 

< ln(r_(-R,R)) > = < ln(r + (-R,R)) > (31) 

and 



< ln(r-(-R, 0)) > = < Zn(r_(0, R)) > (32) 

Consequently, the first two terms in the rhs of Eq.(30) cancel each other and we have only 
to perform averaging of Zn(r_(0, R)). 

These functionals t±(0,R) appear in different physical backgrounds |25|J2"6| . Their 
discrete- x counterpart, which is the sum of products of independent random variables of 
the form 



with 



exp(0 n ), 



is known as the Kesten variable [^] and plays an important role in the theory of renewal 
processes. The distribution function of the continuous-x functional r_(0, R) has been re- 
cently examined in |18| , p8| , p9| , p!9| within the context of diffusion in the presence of a random 
quenched force (the Sinai diffusion f30|,|9|j) and also in the literature on mathematical finance 



The average logarithm of the functional r_(0, R) can be obtained from the probability 
distribution of this functional [28, 29|,19| : 
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2 roc Ah . . 

< ln(T-(0,R)) > = - -r^ (l-exp(-2aRk 2 )7rkcoth{7rk)) - T'(l) - \n(2a) 
i Jo r ^ ' 

/S l -rf?\ 1/2 -| 

« (— ) - r'(l) - ln(2a) + 0(— ), (336) 

where the notation 0(1/ R) means that the neglected terms multiplied by R will give a 
constant as R — > oo. Eq.(33b) has been recently rederived in [[J2|] which tested predictions 
of the replica variational approximation [[J3|] for a particular physical system - a classical 
particle in a one- dimensional box subjected to a random potential which constitutes a Wiener 
process on the coordinate axis |30|j9|| . A detailed discussion of the average logarithm of the 
functional r_(0, R) can be found in f34|| . 

Accordingly, for the DAGS energy we obtain 

E (R) > exp(-4(^) 1 / 2 ), (34) 

7T 

which thus shows that Eq(R) vanishes with an increase of the sample size R not faster than 
a stretched-exponential function exp(— R z ) with z = 1/2. However, this lower bound, which 
is supported by typical realizations of disorder may be improved as we will see in the next 
section. 

IV. LOWER AND UPPER BOUNDS ON E (R). 

In this section we set out to show that, in the limit R — > oo, the dependence of the 
disorder averaged ground state energy E (R) on R is quite different from that in Eq.(34). 
Here we will derive more accurate bounds which show that in the large- R limit, the ac- 
tual dependence of the disorder averaged ground-state energy E (R) on R is described by 
a stretched-exponential function exp(— R z ) but with a smaller exponent, z = 1/3. This 
means that large- R behavior of Eq(R) is supported by atypical realizations of <fi(x). These 
realizations will be also specified below. 
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A. A lower bound. 



Let us begin with the derivation of a lower bound on E (R). We first note that since 
<f>(x) enters the expression for E (R, {0}) only in the form / dx(f)(x), averaging with respect 
to realizations of (p(x) amounts actually to the averaging over different trajectories W(x) of 
a symmetric random walk. Therefore, we can formally write down the average as a product 
of two path-integrals 



where the notations used have the following meaning: The symbol O denotes the set of 
all possible (unrestricted) trajectories W(x) of a symmetric random walk, which "starts" 
at x = at the origin W(0) = 0, and "time" variable x is defined on the interval [0,R]. 
We describe schematically the set f2 in Fig. 2, where for notational convenience we use the 
discrete- a; picture and depict it using the axis W(x) and x, i.e. using "directed polymers"- 
like representation. The trajectories (1) and (2) are two examples of possible trajectories 
which belong to the set f2. The symbol Q' denotes, correspondingly, the set of all possible 
trajectories W'(x) with the "time" variable x defined on the interval [0, —R]. The trajectories 
in O, and Q' are statistically uncorrelated. Finally, the symbols D{W(x)} and D{W'(x)} 
denote that the integration is performed along the trajectories W(x) and W'(x); P[W(x)] 
(or P[W'(x)]) is the corresponding measure of a given trajectory W(x) (or W'(x)), which is 
the standard Wiener measure. 

The next essential step is as follows. Suppose that from the entire set f2 (and Q') we select 
some amount of trajectories having certain prescribed properties and denote this subset of Q 
(fl') as uj (u/). Then, for any positive definite functional E (R, {0}) the following inequality 
holds 



E (R) = < E (R, {0}) > 



[ D{W(x)} [ D{W'(x)}P\W(x)]P\W\x)]E (R,{<l>}), (35) 



f D{W(x)} [ D{W'(x)}P[W(x)} P[W\x)\ E (R,{<j)}) 



> 



> 



[ D{W(x)} [ D{W(x)} P[W{x)\ P[W\x)\ E (R, {</>}), (36) 
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where the integrations in the rhs of Eq.(36) extend only over the trajectories which belong 
to the subsets uj and uj' of the entire sets Q and Q'. Employing the inequality in Eq.(36), 
we get the following bound 

E (R) > [ D{W(x)} [ D{W'{x)} P[W{x)\P[W'{x)\ E (R, {</>}) (37) 

Jul Jul' 

Now we define the subset uj (uj') as follows (Fig. 2): uj (uj 1 ) is the set of all trajectories 
W(x) (W'(x)), which, for any x from the interval [0, R] ([0, —R] for uj'), remain inside the 
strip [—A, A], i.e. such trajectories W(x) (and W'(x)) which obey —A < W(x), W'(x) < A 
for any x from the interval [0, R] ([0, —R] for uj'). In Fig. 2 trajectories which form the subsets 
uj and uj' are exemplified by (2) and (2'). 

Next, we diminish the rhs of Eq.(37), i.e. enhance the inequality in Eq.(37), by substi- 
tuting instead of E (R, {</>}) its minimal value on the subsets uj and uj'. By definition of uj 
and uj', which implies that |W(a?)| < A and |W'(x)| < A we have 

/ dxexp{2W{x)) < 2Rexp{2A), (38) 

J-R 

f° dxexp(-2W(x)) < Rexp(2A), (39.a) 

J-R 



and 

rR 



[ dxexp{-2W{x)) < Rexp{2A) (39.6) 
Jo 



Consequently, for any realization of W(x) or W\x) which belongs to the subs ets uj and uj', 
the following inequality holds 

E (R,{<j>}) > min^{E Q {R,{<t>})} = (40) 

Substituting Eq.(40) into Eq.(37) we find the following bound: 

E (R) > 6XP 2 ( ~ 2 4A) lD{W(x)}£D{W'(x)}P[W(x)}P[W'(x)} (41) 

We notice now that the product of integrals along two "restricted" (statistically independent) 
paths W(x) and W'(x) on the rhs of Eq.(41) is equal to the probability that two independent 
random walkers during "time" R will remain within the strip [—A, A], which means that 
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/ D{W(x)} f D{W'(x)} P[W(x)] P[W'(x)] = P 2 (A,R), (42) 



where P(A, R) is the corresponding probability for a single random walker ]35 



Combining Eq.(41) to (43) we thus obtain 

««> > ^ <| -<-^^>>* m 

The function in the rhs of Eq. (44) contains a free trial parameter A. The inequality in 

Eq.(44) holds for any value of this parameter and thus represents a family of lower bounds. 

Therefore, we will choose such a value of A, which maximizes the rhs of Eq.(44) and thus 

defines the maximal lower bound. For R sufficiently large the maximal contribution to the 

probability distribution in Eq.(43) comes from the term with k = 0, i.e. 

4 n 2 aR 
P(A,R) » -exp(-^£), (45) 



and consequently, 



MR) > J^ 2 exp(-4A - ^) (46) 



Taking the derivative of the rhs of Eq.(46) with respect to the parameter A, we find that 

A = A* = ^(TrVi?) 1 / 3 (47) 

provides its maximal value. Substituting Eq.(47) into Eq.(46) we thus arrive at the following 
"maximal" lower bound 

E (R) > exp(-3(vrV J R) 1 / 3 ), (48) 

which shows that in the large- R limit the DAGS energy vanishes not faster than exp(— R 1 ^ 3 ), 
i.e. at a slower rate that "typical behavior" in Eq.(34). This improved lower bound in 
Eq.(48) is supported by atypical realizations of W(x), such that W(x) oc x 1 ^ 3 , i.e. by 
trajectories of W(x) which are spatially more confined than "typical" realizations of random 
walk trajectories for which W(x) oc x 1//2 . 
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B. An upper bound. 



Let us discuss now the derivation of an upper bound on the DAGS energy. We first note 
that the rhs of Eq.(16) for any given realization of W(x) can be bounded from above 

1 



J-ijrfV)]- 2 dx S« R [^\x>)Y dx> 



+ 



+ lM\x)]-* dx f«M V)] 2 dx> ~ 



1 

+ SM\x)Y* dx ^{Xi)}* dx> 



I-r I-r dx dx' exp(2W(a;') - 2W'(x)) 

-I (49) 

Jo S* dx dx' exp(2W(x') - 2W(x)) 

As one may readily notice, the inequality in Eq.(49) is obtained by simply diminishing the 
limits of integration; in the first term we change the limits of integration over the variable 
x' from [— R, R] to [— R, 0], while in the second one the limits are changed from [—R,R] to 
[0,.R]. Since tp^\x) is positive definite, the diminishing of limits decreases the value of the 
integral and consequently, increases the terms on the rhs of Eq.(49). 

Now we will try to find an appropriate functional of the extremes of the random function 
W(x) which will bound the integrals in Eq.(49) from below, and thus in such a way will 
enhance the bound in Eq.(49). 

We note here parenthetically that this problem turns to be rather non-trivial. In par- 
ticular, standard integral inequalities (such as, for instance, the Schwartz inequality) are 
obviously insufficient since they predict an algebraic growth of the integral 
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r±R r±R 

/ / dx dx' exp(2W(x') - 2W(x)), (50) 
Jo Jo 

while the simple analysis of the "typical" behavior shows that Eq.(50) grows at least as a 
stretched-exponential function of R. In addition, the integrands in Eq.(50) don't possess 
well-defined derivatives and thus one can not expand the integrands in the vicinity of the 
extremes of function W(x) and make use of the standard saddle-point-like estimates. 

To illustrate the derivation of such a bound we first turn to more lucid discrete-space 
picture, assuming that x and y are discrete variables x, y — 0, 1, R, and then approximate 
the integrals in Eq.(50) as products of two sums 

r±R r±R R R 

/ dx dx' exp(2W(x') - 2W(x)) « ^ ^ a xy (51.a) 

-* J ° x=0y=0 

with 

a xy = exp(2W(x) - 2W(y)) (51.6) 

The derivation of the corresponding upper bound in the continuous-space, which is substan- 
tially more lengthy, will be merely outlined in the Appendix to this paper. 

We notice that the rhs of Eqs.(51) is the sequence of (R+l) 2 positive terms, each of which 
is an exponential of the distance between the positions of a given trajectory W(x) taken 
at two different moments of "time" x (summed over all possible x from the interval [0, R]). 
From this sequence of positive terms {a xy } we choose the maximal term, max Xtyt [o,R]{a xy } , 
which is evidently the exponential of the difference of the maximal positive displacement 
M + , (M + > 0), of the trajectory W(x) (which is achieved at some moment x = x*) and 
the maximal negative displacement M_, (M_ < 0), of the same trajectory (achieved at the 
moment y = y*, both x* and y* belonging to the interval [0, R]), 

M + = max xe[m {W{x)} = W(x*) (52) 
M_ = min xe[m {W{x)} = W(y*) (53) 
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Since all a xy > 0, the sum in the rhs of Eq.(51) is evidently larger than the maximal term 
of this sequence, i.e. 

R R 

^2Yj a xy ^ rnax x>ye io tR] {a X y} = exp(2(M + - Af_)) (54) 

x=0y=0 

Eq.(54) represents the (discrete-space) formulation of the desired bound on the integrals in 
Eq.(50). 

Let us now see how this bound can be employed for the derivation of the upper bound 
on the DAGS. Making use of Eqs.(49) and (54) we have that, at a given realization of W(x), 
the ground-state energy can be bounded from above as 

E (R,{<f>}) < exp{-2S') + exp(-2£), (55) 

where we denote by 

S' = max xe[0: _ R] {W\x)} - min xe[0> _ R] {W'(x)} = M' + - M'_ 

and 

S = max xe [ 0iR ]{W(x)} - min xe[0:R] {W(x)} = M + - M_ 

Random variables as S' (or S) are known in the statistics of random walks as a span of 
random walk (see Fig.3), which can be visualized (in d-dimensional space) as the dimensions 
of the smallest box with sides parallel to the coordinate axes that entirely contain the 
trajectory of a random walk |35| . The exact probability distribution P(S,R) of random 



variable S is well-known [pqj ; in the case of large R, a convenient representation reads 



8aR ~ 7T 2 (2k + l) 2 aR , , n 2 (2k + l) 2 aR s , . 

P(S,R) = Et S 2 !] ex P( — — ) ( 56 ) 

Therefore, the calculation of the upper bound on the DAGS energy reduces to the cal- 
culation of the integral 



E Q {R) < 2 dS exp(-25) P(S, R) (57) 
Jo 
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in which, noticing that S and S' have identic distribution functions (although for given 
realizations of W(x) and W'(x) they may have different values) the contribution of 5" may 
be simply accounted by introducing a multiplier 2. 

Let us now consider the behavior of the integral in Eq. (57) in the limit of large R. We 
first note that in this limit in Eq.(56) only the term with k = is relevant. Second, noticing 
that the integrand is a bell-shaped function, we perform the integral using the saddle-point 
approximation. Maximizing the terms in the exponent we get that the saddle-point depends 
on R as 

S = S* = (£^*)i/3 (58) 
2 

and consequently, the bound in Eq.(57) attains the form 

E (R) < 32(^)^ exp( _ _3_ (^1/3) (59) 

Therefore, the upper bound on the DAGS energy shows a stretched-exponential depen- 
dence on R with the characteristic exponent z = 1/3, i.e. aside from the numerical factor 
2^ 1 / 3 in the exponent and pre-exponential multipliers (which are not reliable in view of 
the approximation involved), essentially the same behavior as the lower bound in Eq.(48). 
Since both lower and upper bounds have the same dependence on R and also, on physical 
grounds, E (R) is a monotonically decreasing function of R, we may infer that the stretched- 
exponential dependence with z — 1/3 is the asymptotically exact result for E (R). It is also 
important to note that both the lower and the upper bounds turn out to be supported by 
the same "class" of trajectories W(x), such that W(x) oc x 1//3 . 

To close this subsection we remark that the coincidence in the independence of the 
lower and the upper bounds is, in essence, due to the fact that the measure of the restricted 
trajectories, used in the derivation of the lower bound, and the probability distribution of 



the maximal displacement (or of the span S) are intrinsically related to each other fl3"5 
Actually, the probability P(A, R) that a random walker, starting at the moment R = 
at the origin, remains within an interval [—A, A] in an R-step walk is just the probability 
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that the maximal displacement of this random walker is less than A. It is clear that the 
probability of having the maximal displacement exactly equal to A is given by |35| 

V iA , K) « ^, (60, 

and consequently, the probability of having the span of an i?-step walk equal to S will follow 

P(S , R) „ mm , (61) 

a A \a=s/2 

Further on, evaluating the lower bound we have searched for such an A which maximizes 
the product exp(— 4A)P(A, R). On the other hand, the upper bound was eventually reduced 
to the integral in Eq.(57), which, using Eq.(61) can be rewritten as 

I" dS exp(-2S) P(S,R) w I" dS exp(-2S) dP ^ ^ (62) 

JO JO ' OA \A=S/2 

Integrating Eq.(62) by parts we arrive at performing an integral with the integrand 
exp(— AA)P(A, R). Since the integrand is a bell-shaped function of A and thus the saddle- 
point approximation can be used, calculation of this integral also reduces to maximizing the 
integrand. 

C. Random Gaussian potentials <j){x) with correlations. 

In this subsection we will briefly discuss the behavior of the DAGS in the case where 
fluctuations of 4>{x) are correlated. 

Now, as we have already mentioned the relevant property is the integral W(x) = 
Jq dx' 4>(x'), rather than <p{x) itself. It is therefore convenient to define the correlations 
in random potential 4>(x) in terms of the integral W(x). We consider here the case where 
W(x) is zero in average, as in Eq.(2), and define the second moment as follows 

< W(x)W(x') > ~ \x-x'\ 1+x , — 1 < A < 1 (63) 

The parameter A in Eq.(63) determines the nature of correlations in the random potential 
4>(x). The border-line case A = corresponds to delta-correlated fluctuations of (f)(x), when 



W(x) is a trajectory of the conventional Brownian motion. This case has been examined 
in detail in previous sections. The case of positive A, (A > 0), describes the situations in 
which fluctuations of <f)(x) in two neighboring points x and x' tend to be of the same sign. 
Here the trajectories W(x) have strong persistency; thinking in terms of random walk one 
may say that here the random walker most likely continues the motion in the direction of 
the previous step than changes the direction of motion. Consequently, its trajectories are 
more "swollen" and spatially more extended compared to the case A = 0. Finally, the case 
A < describes disorder with negative correlations when the values of the potential <p(x) in 
two neighboring points x and x' tend to have different signs. Here the random walker has a 
tendency of changing the direction of its motion at each step and its trajectories W(x) are 
essentially more compact, compared to the case of conventional random walk. 

Using Eq.(63) one can readily estimate the typical behavior of the DAGS. Since for the 
typical realizations of W^x) one expects that W(x) ~ x( 1+A )/ 2 , we will have < /nr + (0, R) >~ 
_R( 1+A )/ 2 ; and consequently, 

E (R) typ ~ exp{-R^' 2 ) (64) 

Consider now the behavior of the DAGS stemming from atypical realizations and gen- 
eralize the formalism employed for the derivation of the lower bound. Anticipating the 
reasonings which underly the inequality in Eq.(36) and Eqs.(41), we have that the DAGS 
can be estimated as 

E (R) > exp(-4A) P%(A,R) } (65) 

where P\(A,R) denotes the probability that a random walker, which is at the origin at 
R = and whose trajectories obey Eq.(74) will remain inside the strip [—A, A] during the 



time interval [0, R\. Such a probability can be estimated as p5|j36[] : 

P X {A,R) ~ exp(-R/A d "), (66) 

where = 2/(1 + A) is the "fractal" dimension [36|] of the random walk defined by Eq.(63). 
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Plugging Eq.(66) into the Eq.(65) and maximizing the product with respect to A we obtain 
the following estimate 

E (R) ~ exp(-i? 1 /( 1+ ^)), R > 1, (67) 
or, in terms of the parameter A, 

Eq(R) ~ exp(-R( 1+x ^ 3+x ~>) } R > 1 (68) 

Behavior as in Eqs.(67) and (68) is thus supported by such atypical trajectories W(x) 
which grow with x as a;( 1 + A )A 3 + A ) . It is important to note that again the estimate in Eqs.(67) 
and (68) shows a slower dependence on R as compared to the typical behavior in Eq.(64). 

To close this section we explain what we have in mind when saying that realizations of 
disorder which support the anomalous stretched-exponential behavior of the DAGS share 
common features with the realizations of trajectories which support the anomalous long-time 
decay of the survival probability of a particle diffusing in the presence of randomly placed 
traps or Lifschitz tails in the low energy density of states of an electron in the presence of 
randomly dispersed scatterers [ [T5| , ftSj ]. 

Let us remind, on the example of the trapping problem, some basic formulations and 
results. Suppose a one- dimensional, infinite in both directions, line with immobile traps B 
which are placed completely at random at mean concentration At t = we introduce 
on the line some concentration of particles of another type, say A, and let them diffuse 
independently of each other. As soon as a A particle approaches a B trap, the A particle 
gets annihilated, while the trap is unaffected. The question of interest is to define the time 
evolution of the concentration of A particles (or the survival probability), averaged with 
respect to the spatial arrangement of traps. 

Let C(x,t) denote the local concentration of A particles at the point x at time t. It 
obeys the diffusion equation 

C{x,t) = D ^C{x,t), (69) 
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where D is the diffusion coefficient of A particles. Eq.(69) is to be solved subject to the 
adsorbing boundary conditions imposed at the points occupied by traps; that is 

C(x = X h t) = 0, (70) 

for any Xi from {Aj}, where Xi defines the position of the i- th trap, — oo < % < oo, and 
{Xi} denotes the set of traps' positions. 

A nice feature of the one- dimensional geometry is that this problem can be solved exactly 



15,37|, by simply noticing that evolution of C(x, t) on some interval [Xi, Xi+i] is independent 
of other intervals. Consequently, one has to find the solution of the diffusion equation on a 
finite interval of fixed length W, subjected to the adsorbing boundary conditions at the ends 
of the interval, and then perform averaging with respect to the distribution of the interval's 
length. Such a solution is given by Eq.(43), which in the limit of sufficiently large times 



reads 



Dt 

P(W,t) « exp(-7r 2 — ) (71) 



W 2 . 

Now, the disorder-average concentration of A particles at time t will be defined as 

/•CO 

< C(x,t) > « dW P(W,t) P(W), (72) 
Jo 

where P(W) is the probability of having a trap-free interval of length W. For Poisson 
distribution of traps P(W) behaves as 

P(W) oc exp(-n B W) (73) 

Substituting Eqs.(73), (71) into the Eq.(72) we thus arrive at an integral of essentially the 
same structure as that in Eq.(57), which yields QT3L 



< C{x,t) > w exp ^-3 (^n 2 B Dtf^ (74) 

The behavior as in Eq.(74) shows that the long-time decay of the disorder- average con- 
centration is supported by such bounded realizations W(t) of A particles' random walks 
which obey |W(£)| < A oc t 1//3 , i.e. the same class of trajectories which support the large- R 
behavior of the DAGS in the problem studied in the present paper. 
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V. CONCLUSIONS. 



To conclude, we have studied a new aspect of a one-dimensional localization problem 
associated with the supersymmetric Hamiltonian in Eq.(l) in which the potential <f>(x) is 
Gaussian random function of the spatial variable x. We have derived an explicit expression 
for the ground-state energy of the Hamiltonian (1) defined on a finite interval of the X- 
axis for a given realization of disorder and analysed the dependence of the disorder-average 
ground-state energy on the length R of the interval. We have shown that it is described by 
a stretched-exponential function of the form exp(— R z ), in which the characteristic exponent 
z is dependent merely on the nature of correlations in random potential. In case when 
fluctuations in random potential are delta-correlated we found z = 1/3. In case when 
fluctuations are defined by Eq.(74) we have deduced that z— (l + A)/(3 + A). We have shown 
that such a behavior is quite different from the one expected when only typical realization of 
disorder are considered and thus is supported by atypical realizations of random potential, 
which behave as 



We have also shown that such realizations belong to the class of trajectories which support 
an anomalous long-time behavior of the survival probability of a random walk in the presence 
of randomly placed traps. 



In this appendix, we outline the derivation of the upper-bound for the DAGS in the 
continuous-space limit. Consider the integral in Eq.(50) (for simplicity we suppose that 
limits of the integration are from to +R) and, as it was done before, assume that a given 
trajectory W(x) reaches its maximal value at the point x = x* and its minimal value - 
at the point x = y*. Let us choose some positive constant e, such that < e <C R and 
x* +e, y* + e < R. Since the integrand in Eq.(50) is positive definite, the following inequality 
holds 




Appendix 
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pR pR px* -\-£ py* 

/ dx dy exp(2W(x) -2W(y)) > / dx dy exp(2W (x) - 2W (y)) (AV 

JO JO Jx* Jy* 

Now, taking advantage of the inequality in Eq.(49) we have for the DAGS: 

E (R) < < —5 p > + 

" J R dxJ R dy exp(2W(x)-2W(y)) 



+ < f_ R dxf_ R dy exp(2W(x)-2W(y)) > 
= 2 < f R dxf R dy exp(2W(x)-2W(y)) 

POO p p 

= 2 j dS P(S, R) J J dM + dM^ 6(S - M+ + M_) x 

X < J R dxJ R dyex P (2W(x)-2W(y)) >l<" + =^"-=iW>' ^ 

where the brackets with the subscript (M+ = W(x*); M_ = W(y*)) mean that the average 
is taken with respect to the trajectories W(x) whose maximal positive displacement is equal 
to M + and the maximal negative displacement is equal to M_. 

Further on, the inequality in Eq.(60) enables us to enhance the bound in Eq.(61) and 
write 

Eo{R) < 2 J Q dS p ( S > R ) J J dM + dM - 6 ( S + M_) x 



X K fZ + °dxff +e dye X p(2W(x)-2W(y)) >l^ + =W^-=^)) ^ 
Let us now estimate the value of the following functional: 

E = J^dS P(S, R) J J dM + dM_ 5(S - M+ + M_) x 



X < j::^dxjf^dyex P (2W(x)-2W(y)) >W + =W)*<-=W)) ^ 
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To do this we enclose the points W(x*) and W(y*) by circles of radius 5 (see Fig. 3), where 
5 = S(R) is a slowly growing function. The choice of the dependence S(R) will be made later. 
Further on, we divide the set of all possible trajectories fl into two different subsets. The 
first subset {A} comprises all such trajectories W(x) of random walk (with its maxima at 
M + and minima at M_) which, on the interval x e [x*,x* + e] don't cross the circle around 
the point W(x*) and on the interval x e [y*,y* + e] don't cross the circle around W{y*) 
(e.g. the trajectory 1 in Fig. 3). The subset {B} comprises the rest of the trajectories (for 
instance, the trajectory 2 in Fig.3). We write now 

E = A + B, (A5) 

where A stands for the average of the integrand in Eq.(62) with the trajectories forming the 
subset {^4}, while B denotes the contribution to E coming from the average of the integrand 
over of the trajectories forming the subset B. 

Consider first the contribution coming from the trajectories in the subset {A}. By 
definition of {A}, we have that on the interval [x*,x* + e] the trajectory W(x) obeys the 
inequality M + — 5< W(x) < M + ; and on the interval [y*,y* + e] the trajectory W(y) obeys 
M_ < W(y) < M_ + 5. Consequently, the integrands in Eq.(A4) are bounded from below 
by 

exp(2W(x)) > exp(2M+ -25), (A6) 

exp(-2W(y)) > exp(-2M_ - 25), (A7) 

and thus A is majorized by 

A < eXp( 4 ^ [°°dS P(S,R) exp(-2S) (AS) 
e 2 Jo 

Next we estimate the contribution from the trajectories of the subset {B}. Here, for 
x e [x*,x* + e] and y e [y*,y* + e] the function exp(W^(:r) — W(y)) is always greater than 1 
and consequently 
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r +£ dx ( y +£ dy exp(2W(x) - 2W(y)) > e 2 (A9) 

Jx* Jy* 

Accordingly, the contribution coming from the trajectories of the subset {B} can be ma- 
jorized by 

B < - 2 P({B}), (A10) 



e 2 



where P({B}) denotes the measure of trajectories forming the subset B. When 5 is chosen 
such that 5 2 3> 2 ere, this measure vanishes as 

lnP({B}) ~ (All) 

Now, gathering Eqs.(A8) and (A10) we have that E is bounded from above by 
exp(4<5) - 100 1 x2 



E < 



e 2 



1 \2 

f°°dSP(S,R) exp(-2S) + - exp(-— ) (A12) 
Jo e z 2ae 



Our previous analysis shows that the integral over the span variable S in the first term on the 
rhs of Eq.(Al2) vanishes with R as a stretched-exponential function of the form exp(— R 1 ^). 
Thus, the rhs of Eq.(Al2) behaves as 

- exp(- J RV3) + J_ exp( _iL) m) 

Now we have to make the choice of e and 5(R). One readily notices that the proper choice 
will be if we suppose that e = constant and 5(R) ~ R 1 , where 7 is an arbitrary number 
from the interval ] 1/6, l/3[. If 7 > 1/6 the second term on the rhs of Eq.(A13) vanishes 
with R faster than the first term and thus the leading large-i? behavior will be given by 
the first term on the rhs of Eq.(Al3). On the other hand the requirement 7 < 1/3 insures 
that the leading large- R behavior follows the exp(— R 1 ^) dependence, since R~ 1 ^ 3 5(R) — > 
when R — > 00. 

Therefore, we have shown that also in the continuous-space limit the upper bound on the 
DAGS vanishes with R as a stretched-exponential function with the characteristic exponent 
z = 1/3. The bound derived here (although it suffices to prove the asymptotically exact 
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dependence exp(— i? 1 / 3 )) turns out, however, to be worse than the one found in the discrete- 
space case; it differs from the bound in Eq.(59) by an additional multiplier which grows 
with R as exp(i? 7 ). Besides, this bound is not optimal; there are no well-defined values 
of e and 5 which minimize the upper bound. Apparently, an optimal upper bound in the 
continuous-space can be also devised, but this is beyond the aims of the present paper. 
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Figure Captions 

Fig. la. Schematic picture of a two-level random potential <j>(x). 
Fig. lb. An auxiliary configuration <j>(x). 

Fig. 2. Schematic representation of the sets Q, Q' and the subsets cu, cu'. The set fl (dashed 
triangle on the half-plane x > 0) comprises all possible realizations of an i?-step random 
walk trajectories W(x) with xe[0, R]. The set fl' (dashed triangle on the half-plane x < 0) 
comprises respectively all possible trajectories of an i?-step random walk with xe[0, — R]. 
The subsets u and u' are the areas cut from the sets fl and fl' by the lines W(x) = A and 
W{x) = -A. 

Fig. 3. Maximal positive, maximal negative displacements and the span of the trajectory 
W(x) with x defined on the interval [0, R\. 
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